Temperature Dependence of Zero-Bias Resistances 
of a Single Resistance-Shunted Josephson Junction 
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Zero-bias resistances of a single resistance-shunted Josephson junction are calculated as a function 
of the temperature by means of the path-integral Monte Carlo method in case a charging energy Ec 
is comparable with a Josephson energy Ej. The low-temperature behavior of the zero- bias resistance 
changes around a — Rq/Rs ~ 1, where Rs is a shunt resistance and _Rq = fe/(2e)^ . The temperature 
dependence of the zero-bias resistance shows a power-law-like behavior whose exponent depends on 
Ej I Ec ■ These results are compared with the experiments on resistance-shunted Josephson junctions. 
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Dissipative quantum systems have been studied for 
several decades in various systems 0. Among them, 
resistance-shunted Josephson junctions have been stud- 
ied intensively as a typical system 0, Q- Theo- 
retically, it has been shown that a single resistance- 
shunted Josephson junction exhibits a dissipation- 
induced superconductor-insulator transition at a = 
Rq/Rr = 1, where Rq = /i/4e^ and Rs is a shunt resis- 
tance 0, • This transition is studied for Ej/Ec ^ 1 
and Ej/Ec 3> 1, and is also expected to occur at any 
ratio Ej/Ec, where Ej and Ec are a Josephson en- 
ergy and charging energy of a junction, respectively. Re- 
cently, the phase diagram has been obtained experimen- 
tally for a single resistance-shunted Josephson junction in 
the a — {Ej/Ec) plane ^6,^]. Subsequently, the phase di- 
agram has been obtained also for one-dimensional and 
two-dimensional arrays of resistance-shunted Joseph- 
son junctions. For these higher-dimensional systems, the 
corresponding phase diagrams are modified from the one 
for a single junction. 

In this report, we study the temperature dependence 
of the zero-bias resistance Rq. Analytical calculation 
has shown that Rq behaves as 7^20-2 for a > 1 and 
Ej/Ec ^ 1 [To| . The renormalization group approach 
have indicated that also for any ratio of Ej/Eq, below 
a crossover temperature Tcr, the behavior Rq oc 7^20-2 
may be observed |^. However, there may be deviation 
from this result for Ej ^ Eq above T^r- In order to clar- 
ify this deviation, we study the zero-bias resistance by 
means of the path-integral Monte Carlo method. Similar 
numerical study of a single resistance-shunted Joseph- 
son junction has been performed by Herrero et.al. |l3j . 
They have discussed phase fluctuations along the imagi- 
nary time, clarifying the superconductor-insulator phase 
boundary lies on a = 1 even for Ej ~ Eq- They have 
also obtained the linear resistance on the Matsubara fre- 
quency, i?(iw„). They, however, have not obtained the 
zero-bias resistance on the real frequencies, which can be 



directly related to the experimental data. In this paper, 
we calculate i?(ia;„) first, and obtain the zero-bias resis- 
tance by extrapolation as 



Rq = lim R{iujri 



(1) 



We also compare the numerical results with the experi- 
ments of resistance-shunted Josephson junctions. 

The equation of motion of a single resistance-shunted 
Josephson junction is expressed by a phase difference (f> 
as 



0, 



(2) 



where C is a capacitance, Rg is a shunt resistance, and 
'i'o ^/2e is the flux quantum. The potential energy is 
given as 



y((/.) = -i?jCos0-/|^</., 



(3) 



where / is an external current through the junction. This 
equation of motion is equivalent to a dissipative classical 
motion described by 



• dV 



(4) 



Here, a particle mass M corresponds to fi^/SEc, and 
a damping frequency 7 to l/RgC, where Eq = e^/2C. 
In order to consider the dissipation effect, we use the 
Caldeira-Leggett model , in which, a reservoir consist- 
ing of infinite harmonic oscillators is coupled to the sys- 
tem. In the path-integral formalism, the partition func- 
tion Z is calculated by integrating the degrees of freedom 
of a reservoir as Q 
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V(f){T)e 



(5) 
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where K{t) is a damping kernel given by 



(7) 



Here z/„ = 2-Kn/ph is the Matsubara frequency. 

In order to apply the path-integral Monte Carlo 
method, we introduce a discrete imaginary time Tj = 
h(3j / P, where P is a Trotter number assumed to be even 
and j = 0, 1, 2, • • • , P — 1. Then, the path 0(t) is ap- 



After the 



proximated by a discrete path 0^ 
Fourier transformation =0 + 2 X]r=i ^ cos p , 
6fe sin + ap/2(— !)■', the partition function ||SJ) is ap- 
proximately given as 



Zp = 



5ofF,_P 





p/2-1 

n ' 
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P/2 e-^"«--, (8) 



sin' 5 + «fc ) (flfe + bi 
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where a = Rq/Rs and A is a constant. By this mod- 
ification, the dissipative quantum mechanics is reduced 
to a classical statistical problem. Applying the ordinary 
Monte Carlo method to this effective action, we can sim- 
ulate the dissipative quantum system. Using the linear 
response theory, the zero-bias resistance can be related 
to the correlation function 



dT 



as 



i 



R{iuJn)/RQ = Vl^n 



(10) 



(11) 



where rj — a/27r. The zero-bias resistance is then ob- 
tained by the extrapolation 

Generally, an analytic continuation of the correlation 
function to the real frequencies may have a large noise, 
and is ill-defined. Additionally in the present system, 
the correlation function is not analytic at w = due to 
the presence of a characteristic factor |(jj„| of the ohmic 
damping. Therefore, the Pade approximation, which is 
one of the standard methods of the analytic continuation, 
docs not work well except for moderate temperatures. 
Here, since we are interested only in Riiun = 0), we ex- 
trapolate the numerical data toward = by assuming 
a fitting function. At high temperatures, the numerical 
data of P(i(jj„) can be fitted well to the form 



1 



A\ujJ+B 



(12) 



One example is shown in Fig.n(a'). At low temperatures 
(/3 > 5), the deviation from the form H12|l is observed. 
Then, we fit the data to the form 




1500 2000 



R{iujn)^A' + B'\uj^\ + C'ojl 



(13) 



FIG. 1: Frequency dependence of the Matsubara resistance 
(a) for the high temperature region (/3 = 0.1), and (b) for the 
low temperature region (/3 = 25.0). We take a as 1.5 and the 
ratio Ej/Ec as 1. The insets show the fitting procedure (see 
also the text). 



by using the lowest five frequencies. One example of 
this case is shown in Fig. ^ (b). For moderate temper- 
atures, Rq is consistent with the one obtained by the 
Pade approximation within 10 percent. We expect that 
the method adopted here for the extrapolation is accu- 
rate enough to discuss the temperature dependence of 
Po- 

Before showing the result of the zero-bias resistance 
Po, we study a linear mobility defined as fi\ = /ioPo/Pq, 
where /io — 1/?? is a linear mobility of a free particle. 
The temperature dependence of the linear mobility has 
been studied by Fisher and Zwerger |^ for Ej/Ec <C 1. 
They have showed that for a < 1, as the tempera- 
ture decreases, the scaled linear mobility fJ,\/fJ,o shows 
a reentrant behavior; above a temperature T* it once de- 
creases from I, and increases toward 1 below T*. They 
have also claimed that for a > 1 the linear mobility de- 
creases monotonically as the temperature decreases, and 
becomes zero at the zero temperature. The temperature 
dependence of /J.i/^o in the regime Ej ^ Eq has, how- 
ever, not yet studied. The temperature dependence of 
/ii//io obtained in our simulation is shown in Fig. |21for 
Ej/Ec = 1.0. These results show that the linear mobil- 
ity behaves as predicted by Fisher and Zwerger even for 
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FIG. 2; Temperature dependence of the scaled linear mobility 
for (a) Q — 0.5 and (b) a = 2.0. The ratio Ej/Ec is taken 
as 1. In case of a — 0.5, the linear mobility has a minimum 
value at T = T*. 



Ej ~ Ec- The similar behaviors are observed also for 
other values of Ej / Eq . 

Next, the superconductor-insulator phase transition is 
studied through the temperature dependence of the zero- 
bias resistance i?o- We show Rq as a function of the tem- 
perature for various values of a in Fig. 13 Due to the 
difficulty of the extrapolation, the temperature region in 
which i?o is properly obtained, is limited. Particularly, 
the lack of data at low temperatures makes it difficult to 
distinguish the phase of the ground state. The qualita- 
tive feature, however, can be observed in these figures; 
For a < 1, the result shows the tendency of a reentrant 
behavior, while for a > 1 the result shows a monotonic 
decrease as the temperature decreases. Hence, we expect 
that the phase boundary lies on a = Rq/Rs = 1 even 
for Ej^Ec. 

Finally, we discuss the temperature dependence of Rq 
in the superconducting phase {a > 1). For both a — 1.25 
and a = 1.5, the zero-bias resistances show a power-law- 
like behavior {Rq oc T"^) in the low temperature region 
as shown in Fig. 01 This result disagrees with the an- 
alytic result {Rq oc T^"~^) expected in the low tem- 
perature limit; The exponent u seems to depend also 
on Ej/Ec- This indicates that above a crossover tem- 
perature there exists deviations from the result ex- 
pected in the renormalization group analysis, and that 
Tcr is much smaller than the temperature range studied 



FIG. 3: Zero-bias resistance for various a for (a) Ej/Ec ~ 0.5 
and (b) Ej/Ec ~ 1.0. From the top to the bottom, a is taken 
as 0.5, 0.75, 1.0, 1.25, 1.5 and 2.0. 

in this report. Note that the exponent v may change 
depending on the temperature range, and that v should 
not be thought as a well-defined constant. Here, we use 
the exponent v obtained in the temperature range under 
consideration only for a qualitative comparison with the 
experiments. In Fig. [SJ the exponent v is shown as a 
function of Ej /Ec ■ We observed the tendency that the 
exponent is approximately a linear function of Ej/Eq, 
and that the slope depend on a. 

Let us compare the results obtained here with the 
experiment. Experimentally, the temperature depen- 
dence of the zero-bias resistance Rq has been measured 
to characterize the superconductor-insulator transition; 
In Ref. [i| the system is assigned to an insulator when 
dRo/dT < at low temperatures, and otherwise to 
a superconductor. Systematic analysis of zero-bias re- 
sistance, however, has not been performed for a single 
resistance-shunted Josephson junction 0, Q- In the su- 
perconducting regime (a > 1), the zero-bias resistances 
in these experiments weakly depend on the temperature, 
and do not show a clear power-law behavior. These ex- 
perimental results disagree with the present theoretical 
calculation. The origin of this disagreement is not known. 
We speculate that the environment effect through the 
leads may be important to determine the temperature 
dependence of Rq. Finally, we point out that the zero- 
bias resistance observed in one-dimensional Josephson 
arrays Q obeys the power-law behavior in the supercon- 
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ducting phase. In this experiment, the exponent seems to 
depend Hnearly on Ej/Eq and not on a. Numerical cal- 
culation of one-dimensional resistance-shunted Josephson 
arrays and comparison with the experiment remain for a 
future problem. 



FIG. 4; Power-law-like behavior of the zero-bias resistances 
for various Ej/Ec values at (a) a = 1.25 and (b) 1.5. From 
the top to the bottom: Ej/Ec = 0.5, 0.8, 1.0, 1.2 and 1.5. 
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FIG. 5: Plot of the ratio Ej/Ec vs. the exponent v. 



In summary, we have studied the temperature depen- 
dence of the zero-bias resistance of a single resistance- 
shunted Josephson junction by means of quantum path- 
integral Monte Carlo method. We reproduced the result 
predicted by Fisher and Zwerger even though Ej / Ec ~ 
1. We have also shown superconductor- insulator transi- 
tion around a — Rq/Rs — 1. The temperature depen- 
dence of the zero-bias resistance showed a power-law-like 
behavior in the temperature range under consideration. 
In conflict with the renormalization analysis, the expo- 
nent depends on Ej/Ec- These results have been com- 
pared with the experiment on one-dimensional Josephson 
arrays. 
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